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We study the photon self-energy for the in-medium photon in a three-flavor color superconductor 
in the presence of a magnetic field. At strong magnetic field, the quark dynamics becomes (1 + 1)- 
dimensional and the self-energy tensor only has longitudinal components. In this approximation 
there is no Debye or Meissner screenings at zero temperature, but the electric susceptibility is nonzero 
and highly anisotropic. In the direction transverse to the applied field, the electric susceptibility 
is the same as in vacuum, while in the longitudinal direction it depends on the magnitude of the 
magnetic field. Such a behavior is a realization in cold-dense QCD of the magnetoelectric effect, 
which was first discovered in condensed matter physics. The magnetic permeability remains equal 
to that in vacuum for both transverse and longitudinal components. We discuss the importance 
of the Pauli-Villars regularization to get meaningful physical results in the infrared limit of the 
polarization operator. We also find the covariant form of the polarization operator in the reduced 
(1+1)-D space of the lowest Landau level and proves its transversality. 
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I. INTRODUCTION 

The phase diagram of Quantum Chromodynamics (QCD) has been under intense studies in recent years |l|. At 
extremely high densities and low temperatures color superconductivity (CS) has been established as the favored phase 
[2|. The only known natural environment where this state of matter may occur is the supcrdense core of compact 
stars. The low-temperaturc/high-dcnsity conditions needed for the formation of a CS ground state might be also 
produced in low-energy RHIC experiments or in the planned heavy-ion collision experiments at Nuclotron-bascd Ion 
Collider (NIC A) at JNIR and the Facility for Antiproton and Ion Research (FAIR) at GSI @. 

Both the natural and laboratory circumstances where a CS state could be formed can be accompanied by the 
presence of strong magnetic fields. Compact stars are usually strongly magnetized objects. In the case of magnetars, 
surface magnetic fields can reach values ~ 10 14 — 10 15 G [4|]. Because of the high electric conductivity of the stellar 
medium, during the formation of the neutron star the magnetic flux is conserved. Thus, it is reasonable to expect 
much stronger magnetic fields in the star core where the matter density is much larger. However, the interior magnetic 
fields of neutron stars are not accessible to direct observation, so their magnitudes can only be estimated with the 
help of heuristic methods. Estimates based on macroscopic and microscopic analysis have led to maximum fields in 
the range 10 18 — 10 20 G, depending on the nature of the inner medium, that is, whether it is formed by neutrons 
Q, or quarks jg]. On the other hand, large magnetic fields of the order of m 2 ~ 10 18 — 10 19 G can be produced in 
non-central heavy-ion collisions |3(. If, as claimed in [7j, nontrivial topological configurations can separate different 
charged particles, the existing magnetic field could produce a current that may cause observable consequences in the 
final states. Moreover, it has been shown that a sufficiently strong magnetic field can significantly influence the CS 
state @-|li[]]|0J]. 

At this point, it is worth to recall that despite the original electromagnetic U(l) em symmetry is broken by the 
formation of quark Cooper pairs in the Color-Flavor-Locked (CFL) phase [15| of CS, a residual U{1) symmetry still 
remains. The massless gauge field associated with this symmetry is given by the linear combination of the conventional 
photon and the 8 th gluon fields A^ = cos 6A^ — sin #G 8 [15], LL6[ . The field A^ plays the role of an in- medium or 
rotated electromagnetic field, as the color condensate is neutral with respect to the corresponding rotated charge. 
Thus, a magnetic field associated to A^ can penetrate the color superconductor and not being subject to the Meissner 
effect. 

In the presence of an external magnetic field, a less symmetric realization of the CFL pairing, known as the Magnetic 
CFL (MCFL) phase Q, takes place. The MCFL phase has similarities, but also important differences with the CFL 
phase [8l4l0l ll 21 ] . For example, the ground state has different symmetry and is characterized by two antisymmetric 
gaps A and A^ , instead of just one, as in the regular CFL case Q . The change of the ground state symmetry of 
the three-flavor paired quark matter can be traced back to the rotated charges of the quarks that pair and to their 
interaction with the rotated magnetic field. One should keep in mind that even though all the diquarks are neutral 
with respect to the rotated electromagnetic charges, they can be formed either by pairs of neutral or charged quarks 
with opposite rotated charges. As shown in [8j, the gap A only gets contributions from pairs of neutral quarks, while 
Ab has contributions from both charged and neutral quarks, thus, it can directly feel the background field through 
the minimal coupling of the charged quarks with B. At fields large enough that only the lowest Landau level (LLL) is 
occupied, the field substantially modifies the density of states of the charged quarks and the energy gap Ab becomes 
significantly enhanced by the penetrating field Q. At moderate magnetic fields (henceforth when we say magnetic 
field, we actually mean the rotated magnetic field) the energy gaps exhibit oscillations when eB//i 2 [ly] is varied, 
owed to the de Haas-van Alphcn effect. In the MCFL phase, as well as in the CFL one, the fermion excitations are 
gapped, and the gluon fields acquire masses thanks to the Meissner-Andcrson-Higgs mechanism, but the symmetry 
breaking that gives rise to MCFL leaves a smaller number of Nambu-Goldstone fields, all of which are neutral with 
respect to the rotated electric charge [9|. Hence, the MCFL phase behaves as an insulator, as it has no low-energy 
charged excitations at zero temperature. 

The presence of the magnetic field explicitly breaks the spatial rotational symmetry 0(3) to the subgroup 0{2) of 
rotations about the axis parallel to the field. In a recent paper [12[, we proved that this symmetry reduction has non- 
trivial consequences for the ground state structure of the MCFL superconductor. Performing Fierz transformations 
in a quark system with both Lorentz and rotational 0(3) symmetries explicitly broken, we uncovered new pairing 
channels that favor the formation of new condensates. Of particular interest is an attractive channel that leads to a 
spin-one condensate of Dirac structure Am ~ C757 7 . A gap of this type does not break any symmetry that has 
not been already broken, so it is not forbidden in principle. This condensate corresponds to an average magnetic 
moment of the Cooper pairs in the medium. From a physical point of view, it is natural to expect the formation of 
this extra condensate in the magnetized system because the diquarks formed by oppositely charged quarks will have 
a net magnetic moment that may point parallel or antiparallel to the magnetic field. Diquarks formed by quarks 
lying on any non-zero Landau level can have magnetic moments pointing in both directions, because each quark in 
the pair may have both spins. Hence the contribution of these diquarks to the net magnetic moment should tend to 



cancel out. On the other hand, diquarks from quarks in the zero Landau level can only have one orientation of their 
magnetic moment with respect to the field, because the quarks in the pair have only one possible spin projection. 
This implies that the main contribution to the new condensate should come from the quarks at the zero Landau 
level, an expectation that is consistent with our numerical results |12J |. The new gap was obtained to be negligibly 
small at weak magnetic fields, where the zero Landau level occupation is not significant. On the other hand, at 
strong magnetic fields, the condensate became comparable in magnitude to the original condensates, A and Ab, of 
the MCFL ground state [g|, because the majority of the quarks occupy the zero Landau level in that case. Although 
this new condensate is zero at zero magnetic field, we cannot ignore it even at very small magnetic fields because a 
self-consistent solution of the gap equations with A^O, and Ab 7^ 0, but Am = is not allowed. This is easy to 
understand since there is always some occupation of the zero Landau level, as long as the magnetic field is not zero. 
Thus, once a magnetic field is present, An has to be considered simultaneously with the spin-zero MCFL gaps. The 
Am condensate of the MCFL phase shares a few similarities with the dynamical generation of an anomalous magnetic 
moment recently found in massless QED [17| ■ 

In a recent letter [13| , we reported the peculiar electromagnetic response of the MCFL medium, which exhibits a 
magnetoelectric (ME) effect. The ME effect is the capacity of certain materials to have an electric polarization that 
depends on an applied magnetic field and/or a magnetization that depends on an applied electric field. The ME effect 
has been known in condensed matter for many years (see Ref . [18| for a review) , but to the best of our knowledge it 
had not been previously found in the context of QCD. Neither the CFL [19[ nor the MCFL [l3| phases have Debye 
screening, because no infrared electric screening can be produced by diquark condensates that are neutral. 

In the present paper we provide the details of the calculations that lead to the photon self-energy in the MCFL 
phase with two scalar and one spin-one diquark condensates and use it to obtain the linear electromagnetic response 
of the system in the strong magnetic field region. As it will be shown, in this region the MCFL phase exhibits an 
anisotropic electric polarization, manifested by a different dielectric behavior in the directions parallel and transverse 
to the magnetic field. Thus, the MCFL phase behaves at strong magnetic fields as a highly anisotropic dielectric 
medium. 

The paper is organized as follows. In Sec. II, the expression for the one-loop polarization operator in the presence of 
a rotated magnetic field is obtained in momentum space using Ritus's method. In Sec. Ill, we calculate the full quark 
propagator for the MCFL phase at strong magnetic field. Using this propagator and Pauli-Villars regularization, 
the polarization operator tensor in the strong-field limit is obtained in Sec. IV. In Sec. V, we obtain the electric 
permittivity and magnetic permeability of the MCFL phase from the different components of the polarization tensor 
at strong magnetic field in the infrared limit, po — 0,j>3 — ► 0. In Sec. VI, the covariant structure of the polarization 
operator tensor in momentum space is obtained and the gauge invariance in the strong-field region is verified. Our 
conclusions are given in Sec. VII. Finally, in Appendix A, the transversality of the in-medium photon polarization 
operator in the strong-field approximation is independently proved on general grounds. 



II. THE PHOTON SELF-ENERGY ATB/0 

Let us consider a NJL three-flavor model of massless quarks with Lagrangian density 

C = 4,(1$ + ,i 7 o)^ - G(VT^)(W). (1) 

and baryon chemical potential p. The ground state of this theory is the CFL color superconducting phase. In the 
presence of an external magnetic field the ground state of the theory (fT]) changes into the less symmetric MCFL 
ground state [8| . As mentioned in the Introduction, the explicit breaking of the rotational symmetry by the external 
magnetic field allows for new attractive channels and, as shown in [12(, the MCFL ground state turns out to have 
three independent gaps: two scalar (A and Ab) and one spin-one (Am)- All the diquarks in the MCFL phase are 
neutral with respect to the rotated electromagnctism, so the in-medium electromagnetic gauge symmetry that exits 
in the CFL phase also remains intact in the MCFL one. 

In this paper we are interested in investigating the static electromagnetic linear response of the MCFL supercon- 
ductor in the region of strong magnetic fields. For that, we need to calculate the rotated photon self-energy (i.e. the 
polarization operator for the in-medium (rotated) photon) in the presence of the external rotated magnetic field, B. 
Only the quarks with non zero rotated charge can contribute to the polarization operator. Then, the polarization 
operator for the rotated photon in coordinate space is 



u^( x ,y) 4E Tr [r^fotfro^fo,*)] , (2) 



?=± 



with 



S-( X y)-( G h {X ' V)E (^ y) ) (3) 

the Nambu-Gorkov (NG) full quark propagator of the charged quarks and 

T , = (Qr o \ 

the in-medium electromagnetic vertex. In the above formulas e is the rotated electromagnetic coupling, Q — ± is the 
eigenvalue of the rotated charge operator Q = Qx 1 + lx Tg/y/3 with Q = diag(— 1/3, —1/3, 2/3) for (s, d, u) flavors, 
and T§ = diag(— 1/\/3, — 1/\/3, 2/-\/3) for (b,g,r) colors (for details see [12||). Notice that the polarization operator 
@ gets separate contributions from the positively and negatively charged NG spinors 

•-(&)■ •-(&)■ 

with subindex (±) indicating the sign of the rotated electric charge Q. 
The diagonal and off-diagonal elements in ([3]) are given respectively by 



G %) 



^y)-{K)o^ x ^-nQ) G ho^yKy ^ 

E t^y^- G ho^yK) G t^ x ^ w 

with [G.fi ]~ 1 (x, y) = (ij^dp + eQ^A^ ± p^o)5 i (x — y) the bare inverse propagators for quarks (+) and antiquarks 

(— ), and $ - the gap matrix whose explicit form will be given later. 

We can use the so-called Ritus's method [8], [20, [21[ to transform the polarization operator @ to momentum space. 

In this approach, the transformation is carried out by the eigenfunctions EJ-^ (a;), which correspond to the asymptotic 
states of the charged fcrmions in a uniform magnetic field. Assuming a magnetic field pointing along the positive 

"itten as 

H$+\x) = E p ^(x)A(+)+E l ~ 1 ^(x)A(-), (8) 



X3-direction, the E p (x) can be written as 



E^(x) = El-^-Hx)A(+)+E l p ^(x)A(-), (9) 

Here A(tr) = (1 + iaj 1/ y 2 )/2 are spin up (a = +) and down (<j = — ) projectors that satisfy the following relations 

A(cr)+A(-<r) = l, A(a)A(-a)=0, A 2 (a) = A(a), A(a) 7 f = TfA(a), A(<r) 7 £ = 7^A(-a), (10) 

where || refers to the longitudinal components (0,3) and _L to the transverse ones (1,2) with respect to the direction 
of the magnetic field. The non-negative integer I = 0, 1, 2, ... denotes the Landau levels and the eigenfunctions 

E l p ^\x) = Me- l{poX ° +P2x2+P3x3) D^p^)), (11) 

have normalization constant Mi = (47r|ei?|) 1/ ' 4 /vZT. The functions Di(p(±\) are the parabolic cylinder functions of 



argument pr±\ = y 2\eB\(xi ±p2/eB) and index given by Landau level numbers. 
The E p functions satisfy the equation 

J d 4 xE p (±) (x)E p ', (±) (x) = (2nfm (p p')U {±) (I), (12) 

with E p = 7 o(E p )t 7o , s^(p - p') = S ll '8(p Q - p' )5( P2 - p' 2 )5(p 3 - p> 3 ), and n (±) (Z) = A(±)5 10 + (1 - S 10 ). 



The Ep (±) functions © satisfy 

{id, ± iA^^Kx) = E^Kt^), (13) 

with p( ' given by 



p^ = (po,0 7 ±^2\SB\l 7 p 3 ). (14) 

Since photons are neutral, the transformation to momentum space of the LHS of ^ can be carried out by the usual 
Fourier transform, where we assume translational invariance of the photon self-energy in the presence of a uniform 
magnetic field 



(2ir) 4 5 {4 \p - p'Wip) = / d 4 xd 4 x'e-^ p - x - p '- x ' } U^(x, x'). (15) 

In terms of the E p functions, the charged quark propagator can be written as 

S (Q)^ X ') =$ (^^'(^(giOSlg)^^?^^), (16) 



where 



*<«<*«) - 1 2«tt y^l ) ■ (») 



3 «> < *' G r«'<* ( *) ) 



and 2 ( 27r ?4 = Sz / P (2^) 4?l3 • ^he ma trix ILq-,(?) in (fT6|) ensures that only one spin projection contributes at the 
LLL. 

Hence, in momentum space the rotated photon self-energy reads 



(2.)^) (p - P ')n-b) =f E / ^ 4 *']£ (0^ (0- 



2 —i(p-x-p'-x') 

) 4 



x Tr 



[^E^n^o^*^ 



(18) 



The integrations in the spatial coordinates give [21 1 



/^^«> W r E f>(x,e- = M V»( ? + I >- *-**«->/-.^/» 



and similarly, 



V - T L=e"9 n ^)( n - n '^J„ n -(px)A( CT ) 7 ' i ■A(<t / ), (19) 

— ; Vn!n'! 



/,VEf»(^ h X'*w^-( 2 .)^(, + y-^.'--- e -<- 



x V -i= e ^"(^)( fi - fi '^J rW (p ± )A( ( 7) 7 -A(a'). (20) 



with n = n(l, a), n! = n(l', a'), n = n(l, a) and n' = n{V , a'), defined according to 



n(l,a) = l + sgn(eB)^-^, (21) 



The normalized transverse momentum and the polar angle are defined as p± = yp\ + p 2 and (p = arctan(p2/pi), 



respectively, where p, = p,/ \/2\eB\; the delta function 

5 {3) (q+p- k) = 5(q Q +p Q - k )5(q 2 + p 2 - k 2 )S{q 3 + p 3 - k 3 ), (22) 



and 



min(n,n ) 



Jnn'(P±) 



..In/! 



*-^ m\(n — m)\(n' — m)\ 

m=o v ' v ' 



isgn(eB)p^ 



n-\-n —2m 



(23) 



The presence of the delta functions in (fT9|) and ([20]) enable integrations over ko, &2 and &3 in (fT8|) yielding <5^ 3 '(p— p') = 
<Kpo — Po)^(P2 "^2)^(^3 — P3) an< i the integral over qi gives rise to 5(j>\ —pi), which combined match the delta function 
on the LHS of (0. Tnus ; we have 



e 1 



n""(p)=~|SB|£E£ 



2 

x Tr 



/./' 



j3„ isgn{eB)(n~n' -\-n — n')tp 



e P± Jnn'(p±)Jfin'{p±) 



(2tt) 3 Vn!n'!fi!n'! 

[A(a) 7 ^A(a')n ((5) (/)5^ ) (pW' - ^JA^yA^,^^)' 



(24) 



where [cr] means summing over a, a' , a and a'. Because of the factor e p J- in the integrand, contributions from large 
values of p± are suppressed. This allows us to keep only the terms with the smallest power of p± in J nn '(p±) 21 1 



[max(n,n')]! r .„ 

■/nn'(pj_ "^ — ; 777— M 

97 — 77/1 



Tl-On.n' ■} 



(25) 



Then, the self-energy becomes 



d 3 9 



6 U' kl ^ v 



^ Vn,n' Un,n' 



Q 1,1' W 



x Tr 



[A(a) 7 ^A(a')n ((5) (/)5; (5) (pW' - ^A^A^n^')-^) («<*)] . (26) 



Taking into account that 



5 n ,n' — &l',l5aa' + &l'+a,l&- 



Sn,n' = 8[i,l5a,a' + Sl'+a' ,l8-a> ,a, 



(27) 



we can sum in [a] and /'. Introducing the finite temperature formulation by replacing J dqo — > T J^ , with Matsubara's 
frequencies, <?o = (2r + l)nT, r = 0,±1,±2, ..., we find 



n"»=y|e£|££T£ 

Q l 1° 



dq 3 _ f 

re t 



(2nY 



{Tr[A(+) 7 ^A(+)n (( 3 ) (0)^ ) (F^)A(+)7^A(+)n ((5) (/)^ ) ( g -) 



Tr 



[A(+)y»A(+)n (5) (os^ (F r ff)A(-) 7 "A(-)n (4) (0^ } («) 



Tr 



A(-) 7 ' 1 A(-)n ((5) (05j (5) (F r g) 



x A(+)YA(+)U i0) (l)S[ Q) (q) 



Tr 



A(-) 7 ^A(-)n ((5) (0^ ) (F^)A(-)7 ,y A(-)n (( 3 ) (0^ ) ( g -) 



+ Tr 
x A 



A(+) 7 "A(-)n (c5) (0^ ) (^^)A(-)7^A(+)n (0) (? - l)Sg}(?) 



Tr 



A(-)^A(+)U (C M)S l <d Jp—^) 



L (.Q)v j {Qy 



,(+)YA(-)n {Q) (l + l)S l ± 1 ) (q)]}. 
Here, for simplicity we defined p — q = pv*' — q^> . 



(28) 



III. THE FULL QUARK PROPAGATOR AT STRONG MAGNETIC FIELD 



In order to calculate (|28|). we first need to find the full quark propagator in momentum space S l ,=.(k). Since we 

are interested in the strong-field region, it is enough to find the LLL propagator Sy%?(k), since at strong fields all the 
charged quarks are constrained to the LLL. 



The diagonal 



-l 



G k 1 ^ - fe 1 ]- 1 ® - ^rJmf^ (29) 

and off-diagonal 

s f«)'( s ) = - G T<» , (*)*«) G f<»o , (*)- ( 3 °) 

elements of ffT7|) can be found from the inverse bare propagator [G~= ] _1 (^) and the gap matrix $7^ of the charged 

{Q)0 (Q) 

quarks, with $~^ = 7o[$^ ) ] t 7o- 

The color-flavor structure of the MCFL gap matrix is [12[ 

$?f = Ae^ 3 e lj3 + A B (e^V + e^%- 2 ) + A M [e a/31 (<^ 3 + <^M + e a/J2 (<^ 3 + <M,i)], (31) 



where a, /3 and i,j denote color and flavor indices respectively. Notice that all the coefficients in (|3Tj) arc actually 
matrices in Dirac space, for which we defined A = AC75, A# = A#C75 and Am = AmC75Ci2- In the subspace of 
the positively ipJL) — ( u bi u g, s rd d r c) and negatively ip7_) = i s r, d r ,Ubc> u gc) charged NG quarks, the corresponding 
gap matrices ^t+\ are proportional to the identity in color and flavor and have Dirac structure 

$+ ±) = (-A B ± A M )C 75 A(+) + (-A B t A M )C 75 A(-). (32) 

In the LLL limit (1 = 0), for positively charged quarks we have 

[Gf+v^V 1 ^ 11 ) = (ko ± M)7o - his = 7o Et fc « ~ ( efc 3 T **)]A£,. (33) 



and 



'(+)o 



A 



*r + )[ c r + )o ] $ f + ) - 70 e , n _J: + „^ A * - a ^ 2a (±) + ^ + ^wa^. m 



e=± 



ko - (ek 3 ± /x) 



with A£ = (1 + e7o73k3)/2 being the projectors into states of positive (e = +) or negative (e = — ) energy. At this 
point we concentrate in the propagator for the positively charges quarks, because the calculation is totally analogous 
for the negatively charged quarks. Then, taking into account (|33j) and (|34|) , we can easily obtain G ± and E ± in the 
LLL, given respectively by 

n± '=°/7,!h \r^ fco T (M ~ ek 3 ) ± fc T (^ ^ e k 3 ) , 

G (+) ( fc )=2_. fc 2 _ r e e 12 A k 3 7oA(+) + ^ fc2 _ r-e i2 A k 3 7oAH, (35) 

and 

-f+)"V) = ± E p%p75A^A (+ ) T E ^%p75A^A(-). (36) 

e=± ~ !■ %J e=± ° ~ !■ fc 3 J 



where we defined A = Ab — Am and A = A^ + Am- The quasiparticle energies are el = ^J(p, — ek 3 ) 2 + Aq and 
r k3 = y/& - ek 3 r + A§. 

IV. REGULARIZED POLARIZATION OPERATOR AT STRONG MAGNETIC FIELD 

Now we consider the in-medium photon polarization operator in the strong magnetic field limit, in which all particles 
are restricted to the LLL only. In this limit, the propagators for both positively and negatively charged particles only 
differ in the spin projection of their condensates [L2| ]. Moreover, as we shall see below, the trace over Dirac space 
is independent of the spin projection and thus the contribution to the polarization operator of the propagators for 
positively and negatively charged particles are identical. Therefore, it is enough to calculate the contribution from 
the positively charged particles and multiply it by two. 



In the LLL the charged quark propagators only depend on the Dirac structure gn • 7", with q» = {qo,Q.z) and 



(7°,7 3 ). Taking into account the relation (fTOj) . the polarization operator ((28)) in the LLL reduces to 



IV (p 11 ) = r\eB\T 



go 






Tr 



a(+)7]!^;?(p 11 



gll)A(+) 7 jjS$V) 



(37) 



In (|37|) . the approximation p±^ ~ was already used, since we are mainly interested in the infrared behavior of the 
polarization operator, so the exponential e -23 ^ was eliminated (i.e. e _j5i ~ 1) in (|2"5|) . From (|37p . one can see that 
the theory has been effectively reduced to (I + 1) dimensions. Obviously, the integral (together with the Matsubara 
sum) in (|37|) is ultraviolet divergent by simple power counting analysis. In order to obtain meaningful results, we 
need to regularize Q37p. Because the Lorentz symmetry breaking at finite density, dimensional regularization is not 
appropriate. There is an extra subtlety with using dimensional regularization in this case that concerns to the 75 
dependence of the condensates. As known, 7 s is an intrinsically four-dimensional object with no analog in the higher 
dimensions required in the dimensional regularization approach. Thus, in the following, we will use the Pauli-Villars 
regularization by introducing a counterterm in the polarization operator, which corresponds to free fcrmions with 
zero gap but with a large regulator A playing the role of a mass. At the end, the regulator A will be set to infinity, 
eliminating the introduced degrees of freedom that should not contribute to the physical results. In the Appendix, 
we prove the legitimacy of this regularization scheme in securing the gauge invariance of the polarization operator of 
the rotated photon. Under those prescriptions, the regularized polarization operator reads 



iWp'% 



dqz 



{Tr [At+hJl^V 1 - 9 II )A(+)7JI^(9 11 )] - Tr [a(+) 7 ^a(p 11 - g ")A(+)7^ A (g 11 )] } , (38) 



Here, SaC?") is the propagator of the regularization counterpart, whose inverse reads 



Sa (fc" 



[GtflX 1 ^") 



[G»]l\kW) 



(39) 



with 



[GjlX 1 ^ 1 ) 



&o7o ~ ^37 



A. 



(40) 



Obviously, it is not difficult to calculate the counterpart contribution to the polarization operator. It is actually 
similar to the (real) photon self-energy in QED with the electron mass replaced by A, for which the results can be 
found in the literature [22| • We therefore will not present the detailed calculation of this regularization counterpart 
and only give the final results. 

The trace in (|38[) extends to Nambu-Gorkov, color-flavor and Dirac space. First, performing the trace over Nambu- 
Gorkov space, and using the notation (|I7|) . we have 



n^(pll)=e 2 |e5|T^ 



^'n 



dk 3 
(2tt)= 



rTr 



c,/,s 



1=0 , 



A(+)7l!G+ + ;- U (fc ll )A(+) 7 ||G+ + ;- U (fcll-pll) 



+ At+hJlG-^-V^At+hllG^-V 1 -P^ 

+ A(+)7|!s+ + ; =0 (fcll)A(+)7j|S- + / = °(fc l1 -p«)] +LV(A). 



P"( + ) (^1) 



Where LT^ (A) corresponds to the contribution of the second term in the RHS of ([38]) . 

Let's introduce the mixed representation for the quark propagator in the LLL (sec Ref. |23| for a similar represen- 
tation at zero field) 



G 



(+) 



V,fc3)^T£e- fe -Gf + /= U (fcl'), 



(42) 



and its inverse 



G 



(-) 



{*$)=[' dre fc ^Gf + /-°(r,fc 3 ). 



(43) 



Then, after performing the Matsubara sum in terms of a contour integral in the complex fco plane, we have 

G+ + / = V k 3 ) = £ A£ a7o {K 3 [6{-t) N(ei 3 )] e e V (1 - < 3 ) [6(r) - N (<£,)] e"<. T } , (44) 



G- + /- (r, fc 3 ) - - $> A k3 {< [0(r) - JVC4,)] e- e V - (1 - <J [0(-r) - AT(e k3 )] e*. r } 



(45) 



with JV(x) = (e x ' T + 1) _1 and n k = (e k + /i — efc 3 )/2e£ being respectively the fermion distribution function and 
occupation numbers of particles (e = +1) or antiparticles (e = —1) at zero temperature. Substituting (|42[) into (|41|) . 
and making use of the Kubo-Martin-Schwingcr relation for fermions 



Gf + /- (l/T-r,k 3 ) = -Gf + ^'(-T,k 3 ) 



y± 1=0, 

'( + ) 



and the identity 



T^ koT = E (-ir^-S 



?n— — oo 



(46) 



(47) 



we obtain 



T^> s A(+) 7 HGf + ) (fcll)A(+) 7 IJGf +) (g 



l/T 



dre Por T 



A(+)7l!Gf + /-°(r,fc 3 ) 



x A(+)7jJ7oGf +) (r,g 3 )7o 



(48) 



Here, we used that e p °/ T = 1 for bosonic Matsubara frequencies pa = —i2mT and defined q^ = k^—p^. Now inserting 
(|4"4" |) - (|45p into (|4"5)l and integrating over r, we have 



T^Tr s [A(+) 7 llG+ + ;- (fcll)A(+) 7 ,!G+ + ;- ,, ( (Z ll 

fco 






Po + e£ 3 + e«; p - e^ 3 - e q ' 3 



l-i\T(e? )-7V(e q J 



(l-n^)(l-^ 3 ) njn^ 



VQ - 4 3 + e q 3 Po + 4 3 - € 



Q:l 



*(£,)"*(£,) 



(49) 



and 



r^Tr s [A(+) 7 llG (+ ;-°(fcll)A(+) 7 ,!G ( ;;- U ( g ll 

fco 

' (1 - nt )n e ' nt (1 -n e ')\ 

v K3 1 q 3 k 3 V q 3 / I 



Ev, 



l P0 + 4» +e q 3 P°- £ k 3 - e q 3 



l-JV(e k3 )-7V(e q J 



"k 3 n q 3 ( 1 -"k 3 )( 1 -"q 3 ) 



\Po - % + e q3 po + ^ 3 - e q3 

In (l49l)-(l5"0l) we introduced the notation 



N(eU-N(ei 3 ) 



T* — TV 



7 o 7 |lA± e A(+) 7 o 7 l!A± e 'A(+) 



(50) 



(51) 



The same procedure can be applied to the off-diagonal components in (|41[) to end up with 

T^Tr[A( + ) 7 ^ +) ^(fcll)A(+) 7 J|sf +) ^ ( g ll); 



= -£*£ 



i 



4e k 3 £ q3 1 V P0 + £ k3 + <3 PO ~ £ k3 - < 3 



l-iV(eL)-iV^ 3 ) 



Po-4 3 +< 3 Po + e k3 -e q ' 



(52) 
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with 



U% = Tr s 



7 5 A± e A(+) 7 l| 7 5A^'A(- 



(53) 



Putting everything together, the polarization operator flU} for rotated photons reads 

dk 3 



iW) 



2£ 2 \£B\ 



(2n)'- 



E 



ni(l-n 2 ) (1 -ni)n 2 ^ ^ _ ^ _ ^ , + /'(l -ni)(l - n 2 ) 



Po + £i + £2 Po - d - «2 
(l-ni)n 2 ni(l-n 2 ) 



nin 2 



Po + (i + e 2 Po - ei - £2 
1 1 



Po + (i - £2 Po - ei - £2 

A 2 

2 +IV A 

4eie 2 



(1 - Ni - N 2 ) + 



(1 - Ni - N 2 ) - 



PO - £l + £ 2 PO + £l - £2 

mn 2 (l-m)(l-n 2 ) 



Po - £l + £2 

1 



PO + £l - £ 2 

1 



{Ni - N 2 ) 
(JVi - iV 2 ) 



T" 1 

' //J 



' /if 



PO - £l + £ 2 Po + £l - £2 



(JVi - 7V 2 ) 



x (W+ +W-) 



(54) 



Here some compact notations, similar to those used in [23[, were introduced: 1 stands for subscript k3 and superscript 
e, and 2 for subscript q3 and superscript e'. Notice that the prefactor 2 comes from the trace in color-flavor, since 
the RHS of (|4Tj) is diagonal in color-flavor space. 

The remaining trace over Dirac space in (|51[) and (|53l) can be easily performed to obtain 



77, 



± 
00 



-^00= 2 (1 + e6 ' k3 ^ 3) ' 



(55a) 



% 



± 

03 



'30 



w 03 = W 30 = o(=Fek 3 T e'q 3 ), 



(55b) 



7^=Z4 = -(l + ee'k3q3). 



(55c) 



Notice that we would end up with the same results if we replace the spin projectors A(+) in (fSTj) and ([53| by A(— ). 
That shows that the contribution of the positive and negative charged particles with corresponding spin up (A(+)) 
and spin down (A(— )) projectors respectively, to the polarization operator are identical, as we discussed previously. 
The different components of the rotated-photon polarization operator (|54j) can be explicitly given as 



iwp"; 




po - £1 + £2 
x {l-N El -N E2 ) 



po + £l - £2 
E 1 E 2 -k 3 g 3 -A 2 | 

E1E2 



( - 

Po + £1 + £2 Po 



£1 — £2 



-Xl-Nt-Ni) 



£l£ 2 - gi&z - Ag 
2£l£ 2 



2£l£2 
1 



1 



po + E x + £ 2 
1 



po — E\ 4- £7 2 po + -Ei - £"2 



Po 

(N El -N E2 ) 



Si — B 2 
^i^ + fcgga + A 2 



E\E 2 



(56) 



iWp") 



e 2 |eB| 



<Jfc3 



^(ek 3 + e'q 3 ) 



PO + £1 + £ 2 po - £l - £2 



^-k 



— ^— +^ — W-Wf+f- 

p -£l+£ 2 P0 + £l-£2/ \ 2e l 2e 2 



x (l-JV Bl -Ne 2 ) 



<!■•> 



h 



E 2 Ei 



1 



+ l 

1 



1 



1 



po - E x + E 2 V q+E x -E 2 ) \E 2 Ei 



po + Ei + E 2 po - Ei - E 2 
Q3 + ^) (N El - N E2 ) 



(57) 



and 
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iWp") 



e 2 \eB\ 



dk 3 



1 



^(l + ee'kaQs) 



1 



E X E 2 - k 3 q 3 + A 2 



x {l-N El -N E2 y- 



E1E2 



(JVi - JV 2 ) 



Po + ei + e 2 Po - ei - £2 



eie 2 +£i6 -A§ 



2ei£ 2 
1 



1 



(1 - JVi - iV 2 ) 



£ie 2 - 66 + A 2 , 



2eie 2 



1 



1 



Po + E r + E 2 po-Ei- E 2 

EiE 2 + k 3 q 3 - A 2 



Po- Ei + E 2 po + Ei - E 2 



(N El -N E2 y- 



E\Ej2 



(58) 



Here, we used the notation E\ = \Jk\ + A 2 , E 2 = y/ ' q\ + A 2 and 6 = ek 3 — /x, 6 = e'93 — M' 



ELECTRIC PERMITTIVITY AND MAGNETIC PERMEABILITY IN THE STRONG-FIELD 

APPROXIMATION 



In the following, we will consider the zero temperature limit, T = 0, of the polarization operator components 
(|56| - (|58|) . and consequently all the Fermion distribution functions vanish. Moreover, we can drop the contribution 
of the antiparticles because it is always much smaller than that of the particles. Because we are interested in the 
infrared behavior of the photon self-energy, we expand (|56[) - (|58p in powers of the photon momentum components po 
and p 3 , up to quadratic terms 



n, 



00 B 



= — lim 

A— KX 



gjeglg| 
6tt 2 



1 

A 2 



J_. 

A 2 ' 



c 2 



eB\p% 



6^ 2 A 2 



(59) 



22| S 



Il33fl = - lim 

A— >oo 



eB\pl ( 1 



6tt 2 



( A^~ 



J_. 

A 2 ' 



~e 2 \eB\P 2 o 
6^ 2 A 2 ' 



(60) 



and n 30 ^ = n 3^ ~ 0. As we expected, the regulator introduced through the Pauli-Villars regularization scheme 
does not appear in the final results once we take A — > 00. 

Because noo has no constant contribution in the infrared limit po = 0>P3 —* 0, one immediately concludes that there 
is no Dcbye screening in the strong-field region, as it was the case at zero field in the CFL phase [l{|. This is simply 
because all quarks are bound within the rotated-charge neutral condensates. There is also no Meissner screening 
(i.e. 033 is zero in the zero-momentum limit), as it should be expected from the remnant U(l) gauge symmetry. 
However, the condensates have electric dipole moments and could align themselves in an electric field. Hence, this 
should modify the dielectric constant of the medium. Since the quadratic term in the effective U(l) Lagrangian is 
given by A^ (— p) [-D^ 1 (p) + n M „(p)]A„(p), with D _1 being the bare rotated photon propagator, the effective action of 
the £/(!) field in the strong-field region is given by 



S, 



eff 



dM-lE|| 



Eii 



-E 1 ■ E 1 



2An 



-H„ H 



1 
2\\ 



H H 



(61) 



where the separation between transverse and longitudinal parts is due to the 0(3) —> 0(2) symmetry breaking 
produced by the strong magnetic field B. In (|61|) . E, H arc weak electric and magnetic field probes, respectively. In 
(|6"Tj) the coefficients e and A denote the electric permittivity and magnetic permeability of the medium respectively. 

From (|5T))) - (f6T)|) it is straightforward that in the infrared limit the transverse and longitudinal components of the 
electric permittivity and magnetic permeability become 



S2| 2 



A 



A 



1, 



!•: e ll = 1 + XMCFL = l 



eB\ 



6ir 2 A 2 ' 



(62) 



where Xmcfl ^ s ^ ne longitudinal electric susceptibility. Notice that the longitudinal electric susceptibility is much 
larger than one because in the strong- magnetic- field limit eB ^ Ag (T3| ■ 

Although a static £/(l) charge cannot be completely Dcbye screened by the U(l) neutral Cooper pairs, it can still be 
partially screened along the magnetic field direction because the medium is highly polarizablc on that direction. This 
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is due to the existence of Cooper pairs with opposite rotated charges Q that behave as electric dipoles with respect 
to the rotated electromagnetism of the MCFL phase. Moreover, the electric susceptibility depends on the magnetic 
field. When the magnetic field increases in the strong-field region, the susceptibility becomes smaller, because the 

coherence length £ ~ 1/A decreases (i.e. A increases) with the field at a quicker rate than y eB pjj, and the 
pair's coherence length £ plays the role of the dipole length. Hence, with increasing magnetic field the polarization 
effects weaken in the strong-field region. The tuning of the electric polarization by a magnetic field is what is called 
in condensed matter physics the magnetoelectric effect. We have discussed this phenomenon in highly magnetized CS 
in detail in [13| . From (|62p . we also see that at strong magnetic fields the medium turns out to be very anisotropic. 
The fact that the electric permittivity is only modified in the longitudinal direction is due to the confinement of the 
quarks to the LLL at high enough fields. 



VI. COVARIANT STRUCTURE AND GAUGE INVARIANCE OF THE POLARIZATION TENSOR 



The photon polarization operator should be gauge invariant. That is, in the strong-field approximation, it should 
satisfy the transversality condition in the reduced (1 + 1)-D space (pjJ'nju,(pH) = 0). As known, the polarization 
operator tensor can be expanded in a superposition of independent transverse Lorentz tensors. The number of these 
basic transverse tensors depends on the symmetries of the system under consideration. For example, in vacuum, where 
the only available tensorial structures are the four-momentum and the metric tensor, there is only one gauge invariant 
structure. When a medium is under consideration (i.e. at finite temperature or finite density), since the Lorentz 
symmetry is broken, there is an additional gauge invariant structure that can be formed by taking into account a 
new four- vector, the four- velocity of the medium center of mass, m m , [24| . When a magnetic field is applied on that 
medium, then the structure of the polarization operator is enriched by an additional tensor, F),„. Then, at finite 
density and in the presence of a magnetic field, there are nine independent gauge-invariant tensorial structures [25j. 
At strong magnetic field, when the particles are confined to the LLL, due to the fact that the transverse momentum 
is zero, there is a dimensional reduction leaving only the tensors g]ii/,pj! and uji = (1, 0) at our disposal. The original 
nine structures of Ref. [25| now reduce to only two 






(P ll fg, 



fa/ 



PM 



and 



T (2) 



II _ pU u1] -p") 

'" (pll) 2 



jl _ pl{v}\ -pH) 



T r 



(63) 



(64) 



Moreover, one can readily check that the two tensors (|63|) and (|64|) are equivalent, which indicates that the rotated- 
photon polarization operator tensor, at strong magnetic field, only has one independent structure 



njUp") = n( P 'i, /.,£?) {p*Y9%-pIvI 



(65) 



with n(p",/i, B) being a scalar coefficient depending on the photon longitudinal momentum, baryonic chemical po- 
tential and magnetic field. 

Contracting with the tensor g^ both sides of Eq. (|65|) . we will obtain the scalar function II (p" , /z, B) by substituting 
Eqs. (|3"5"|) - ([3l)|) into Eq. (|3"?| and performing the traces in all spaces, 



n(piU,B)=(<r)Unii 



fJLV 



4g2 |^E/7^E(l + -'q 3 k ; 



)x 



A 2 

-V) 



Then, after performing the Matsubara sum, we end up with 
n( p ll, M ,3)=-4eW^£(l + 

e,e ; 

l l 



[5o 2 -fe] 2 )(fco 2 -[<] 2 



(66) 



A 2 



4e^ tt 

13 K 3 



PO 



'/;: 



1 



PO 



•</:i 



1 



.^o + eL-e? Po -e°-el 3 



x(l-N(et)-N(eU) 



x(N(e*)-N(el)) 



(67) 
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Now, if we expand the components of the polarization operator (|65|) . after substituting the coefficient H(p, /j,, B) with 
(|57]). to quadratic orders in powers of po and ^3, we will obtain exactly the same results as those of Eqs. (|59p and 
(|60|) . Notice that we do not have to regularize the scalar function (|67[) . because the divergent terms cancel out when 
taking the trace in (|66p. Hence, the medium electric permittivity and magnetic permeability obtained from (1591) - (I00 |) . 
have a gauge-invariant nature, as they are obtained from the gauge covariant form of the polarization operator tensor 
d65l). 



VII. CONCLUDING REMARKS 

In this paper we investigate the polarization effects that affect the in-medium photons in the strongly magnetized 
MCFL phase of CS. With this purpose, we studied the self-energy of the rotated photons associated with the unbroken 
U(l) symmetry of the MCFL phase in the presence of a very strong magnetic field. In calculating the one- loop 
photon polarization operator we used the Pauli-Villars rcgularization scheme to regularize the diagram's ultraviolet 
divergencies. The use of this regularization scheme was crucial to get rid of some unphysical results. We call attention 
that, without such a convenient regularization we would have ended up with a constant Meissner mass, which is of 
course in contradiction with the remaining £/(l) gauge invariance of the theory. This simply indicates that it is not 
allowed to take the infrared limit of the Feynman's diagram while it is not regularized. The regularization procedure 
in this case is however more subtle than in vacuum; where it is possible to use dimensional regularization, hence 
preserving the gauge symmetry of the polarization operator. In CS, nevertheless, dimensional regularization is not a 
suitable regularization procedure, since in this high-dense medium, Lorentz-symmetry is broken, and besides, there 
is an extra complication due to the presence of 75 in the Cooper-pair condensates. As known, 75 is an intrinsically 
four dimensional object that cannot be generalized to higher dimensions. We thus use the Pauli-Villars regularization 
scheme. This regularization procedure is not in general gauge invariant for diagrams with internal lines of gauge 
bosons, since the addition of an auxiliary boson propagator with a mass term (i.e. given by A in (|40p) will break the 
gauge symmetry. However, for CS in the hard-dense-loop approximation, the diagrams are formed only by internal 
fermion lines, and then the Pauli-Villars regularization becomes an appropriate gauge invariant approach. 

Notwithstanding the difference between the photon polarization in magnetized QED in vacuum and the one we are 
reporting for the color superconducting medium, we found that there is still one similarity regarding the covariant 
tensorial structure of the polarization operator in both cases. In the QED case [22}, the polarization operator in the 
strong-field limit exhibits only one transverse structure in the subspace of the longitudinal momenta. In the MCFL 
phase of CS, although for a magnetized dense medium there exist in general nine possible independent covariant 
structures J25J, in the strong-field limit they reduce to only one in (1 + 1)-D, which is similar to that of QED (see Eq. 
(|65|) ). The fact that we can find a covariant transverse representation for the in-medium photon polarization operator 
is another manifestation of the gauge invariance of the results we are reporting. 

We found that there are no Debye and Meissner screenings in the strongly magnetized color superconducting 
medium. The reason for the lack of Debye screening is simply because all the particles in the medium are condensed 
in Cooper pairs that are neutral with respect to the rotated charge, and thus, cannot screen a static rotated charge 
in the medium. On the other hand, since the rotated U(l) gauge symmetry remains unbroken, we should not expect 
a different from zero Meissner mass. Nevertheless, the medium has a large electric susceptibility, which indicates 
that the medium is highly polarizable. A similar conclusion was found for the gluodynamics of the remnant SU(2) 
symmetry in the 2SC phase of CS at zero magnetic field in Rcf. [26|, and for the rotated electric polarization in 
the CFL phase at zero magnetic field [l{|. Nevertheless, the electric susceptibility we found in the MCFL phase is 
different from those phases in two aspects. On one hand, it can change with the applied magnetic field, and on the 
other hand, it is anisotropic. The electric susceptibility in the strongly magnetized MCFL is quantitatively smaller 
than that of the CFL phase .19] for the same value of the baryonic chemical potential. It is due to the fact that the 
electric susceptibility is proportional to |ei?|^ 2 , and the coherent length, £, decreases as 1/Ao, and Ao increases with 
the field at a quicker rate than \JeB. Hence, as the coherent length of the Cooper pair plays the role of the electric 
dipole moment length, when the magnetic field increases in the strong field region, the susceptibility becomes smaller. 
The fact that the electric susceptibility depends on the magnetic field is an effect analogous to the magnetoelcctric 
effect known in condensed matter physics. In [l3j we discussed in details the realization of this magnetoelcctric effect 
in CS, as well as its possible implications for astrophysics. 

Comparing our result for the electric susceptibility (|62p with that obtained in strongly magnetized QED, Xqed = 
{a/3>ir){\eB\/ra 2 ) [22], we find that the two results are quite similar after the following replacement of the dipole 
length 1/Ao — ¥ 1/ttt,. Also, there is an extra factor 2 in X M cfL' wn ich comes from the contribution of two electric 
dipoles associated to the Cooper pairs formed by the charged quarks present in the MCFL phase: (s r , Ub) and {d r , u g ). 
Nevertheless, there is a substantial difference between the QED electric susceptibility and that of the MCFL phase 
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of CS regarding their dependence on a magnetic field. While in the QED case, Xqed increases with B, since the 

coherence lenth, £,qed ~ l/ m i does not change with B, in the MCFL phase, Xmcfl decreases with B, as it was 
discussed before. It is also interesting to see that the chromo-electric susceptibility in strongly magnetized QCD, 
Xqcd = ( a s/6Tt) ~^2 q ^i{\e.B\/m 2 ) [27| . also has the same qualitative behavior as Xqed> with a coefheient 6 instead of 
the 3 found in QED, which is related with the trace in color associated to the gluon vortices. 

The anisotropic nature of the electric susceptibility in the strongly magnetized MCFL phase, together with the 
lack of Debye screening, implies that electric fields transverse to the applied magnetic field are not modified at all 
in the color superconducting medium. Possible effects of this anisotropic susceptibility for compact astrophysical 
objects were discussed in Ref. [13[ . Also, the interplay between the electric properties of dense media and strong 
magnetic fields could be of interest for future low-temperature/high-dcnsity heavy-ion collision experiments where 
high magnetic and electric fields can be generated [3| and presumably CS can be realized [ljl UM ■ 

It is important to point out that possibly the situation will be slightly different in the 2SC superconductor in a 
magnetic field. Even though we expect a similar anisotropy to be present in the 2SC case too, the charged blue 
quarks can in principle Debye screens there an external electric field in all directions. It will be worth investigating 
this case in more detail since 2SC superconductivity can be the more reliable phase at moderate densities and low 
temperatures j29j . 
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Appendix A: Transversality of the rotated photon polarization tensor in the strong-held Approximation 

The rotated electrodynamics in a color superconducting medium is a U(l) gauge theory. As a consequence, the 
in-medium photon self-energy operator should be transverse with respect to photon four-momentum. In the strong- 
field approximation, when all the quarks are confined to the LLL, the photon self-energy tensor is reduced to the 
(0,1) plane of the longitudinal-momentum (1+1)-D space, as has been shown earlier. Although from the covariant 
structure of LL^,, (|65p . the transversality of the in- medium photon self energy is already guaranteed, we want to show 
in this Appendix that the transversality is intimately associated in this case to the structure of the quark propagator 
in the LLL. Also, we want to highlight the relevance of the Pauli-Villars regularization in securing the transversality. 

In the strong-field approximation, where the quarks are confined to the LLL, the inverse quark propagator is given 
by 

[«5(+)°]- 1 (fc 11 ) =5- 1 (fc«) = j&" +M7o^3 + i75A A(+)a 2 , (Al) 

where the Pauli matrices, ctj, act on the Nambu-Gorkov space. With the same compact notation, the vertex (Q| for 
the positively charged particles is given by the diagonal matrix in Nambu-Gorkov space 

T% = rf5 AB - (A2) 

Where A, B are Nambu-Gorkov indices. From (|A1I) we have 

^=S- 1 Qc^+q^)-S- 1 (k^), (A3) 

Sandwiching (|A3|) between <S(fc" + qW) and S(k"), we find that 

S(fc" + q ll )^S(k^) = S(fc") - S(0 + g 11 ). (A4) 

The transversality of the rotated photon self-energy tensor requires (q^y'H^, which from (|37|) is equivalent to 

T E / T^ Tr f 5 ( fc " + <?" )A(+)S0l'S(fcll ) A(+)e 7 i!l = 0, (A5) 

■70 J [Z7r > 

From the fact that [A(+),S(fcH)] = 0, it follows from (|A4|) and (|A5|) that 

r E/(0i{ Tr [ A(+)5(A: " )p: ^] -Tr[A(+)5(fcll+ g ll)g 7 Jl]}-0. (A6) 
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Naively, this requirement is satisfied by simply shifting the momentum. However, the shift of the integration variable 
is no longer legitimate if there is a linear or higher divergence. We thus introduce the Pauli-Villars regularization and 
the regularized version of the LHS of (|A5j) reads 

dq3 
(2tt)S 



T E / (S2 { Tr [ 5 ( fc " + g ll )A(+)^ ll 5(fc l! )A(+)g 7 l|] - Tr [<S A (fcH + gll)A(+)e^ ( S A (fcll)A(+)e 7 l|] } . (A7) 



where S\(k) is the Pauli-Villars propagator with zero gap, zero chemical potential and mass A. We observe that 
5a (fc) satisfies the same relation (|A4|) . Therefore 



c 2 



(ef)llnjL fl = -(e5)T£ 



dq-3 



F v (k^-F v (kW+qW) 



(A8) 



where 



F v (k) =Tr A(+)5(fc")e T 



2JI 



-Tr 



A(+)S A (0)ejS, 



(A9) 



It can be shown explicitly that for a fixed cutoff A, neither of the integrals in (|A8[) is more divergent than logarith- 
mically. Then, the shift of the integration variable is legitimate and we thus obtain 



(g")Hn 



!xv R 



0. 



(A10) 



as required by gauge invariance. 
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